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ABSTRACT 

Infraied  Imagery  from  polar-orbiting  satellites 
ptovides  a  synoptic  and  long-term  view  of  lead 
patterns  in  Arctic  pack  ice.  The  large  quantity  of 
satellite  data  in  image  form  suggests  the  use  of 
Automated  methods  for  compiling  lead  statistics  from 
imagery.  A  Hough  transform  technique  for  the  semi- 
automated  extraction  of  lead  orientation  and  spacing 
is  described  and  first  results  are  presented. 
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1.  INTRODUCTION 

Sonsnts  cn  po la r - or b i t i ng  satellites  produce  digital 
Imagpty  vlth  which  the  regional,  seasonal,  and  annual 
variability  of  large-scale  lead  patterns  in  Arctic 
pack  ice  can  be  studied.  The  infrared  AVHRR  sensor 
on  boat  d  the  NOAA  series  is  well  suited  to  this  task 
because  of  its  nearly  arctic-wide  swath  and  revisit 
frequency  of  several  times  per  day.  Thermal  data 
offns  day  and  night,  year-round  coverage.  However, 
cJnuH  contamination  is  often  a  problem  at  high 
latitudes,  and  because  AVHRR  imagery  has  a  resolution 
of  1km  at  nadir,  only  large  leads  are  resolved. 

Within  these  limitations,  characteristics  of  a  lead 
ensemble  such  as  average  lead  width,  length,  spacing 
and  orientation  can  be  obtained  from  imagery  and 
related  to  atmospheric  and  oceanic  stress  fields  and 
to  surface  fluxes. 

Satellite  remote  sensing  provides  the  spatial  and 
temporal  coverage  necessary  to  compile  a  climatology 
of  l*>ad  character  1  st  ics.  With  this  advantage  comes 
the  task  of  extracting  lead  features  from  large 
amounts  of  image  data.  Here  a  method  for  the  semi¬ 
automatic  extraction  of  lead  orientation  and  spacing 
from  binary  images  of  leads  is  presented.  Although 
demonstrated  with  AVHRR  imagery  the  technique  is 
suitable  for  any  digital  Imagery  in  which  lead  pixels 
can  be  distinguished  from  background  pixels. 

2.  THE  HOUGH  TRANSFORM 

The  Hough  transform  is  a  computer  vision  technique 
for  detecting  lines,  circles,  or  other  shapes  in 
imagery.  Recently  it  has  found  application  in  remote 
sensing  (Cross,  1908),  where  its  speed,  flexibility 
in  extiactlng  shape*;,  and  relative  insensitivity  to 
Image  noise  makes  it  an  attractive  alternative  to 
other  methods  for  object  recognition.  To  implement 
the  ttnnoform,  an  edge  detector  or  other  method  for 
identifying  leads  Is  first  applied  to  the  image. 


Each  pixel  of  the  image  is  classified  as  either  a 
lead  pixel  or  not  a  lead  pixel.  Fixels  identified  as 
lead  pixels  in  image  space  are  then  mapped  into 
"parameter  space",  where  parameters  are  those  vhi«*h 
describe  the  shape  being  sought.  Here,  it  Is  assumed 
that  lead  pixels  are  roughly  collnear  points, 
therefore  the  transform  is  employed  as  a  line  finder 
using  the  normal  parameterization  of  a  line 
p  =  x  cos0  +  y  sin0 

where  0  is  the  angle  of  a  normal  to  the  line  and  is 
its  distance  from  the  origin  (Fig.  la).  Farameter 
space  is  represented  as  an  accumulator  array  of 
discrete  O.p  values.  For  an  (x,y)  image  pixel,  9  is 
incremented  from  0  to  100  1 n  steps  of  one  degree  and P 
calculated  for  each  9  .  Accumulator  element  (6.p)  is 
incremented  for  each  pair.  Every  lead  point  In  the 
image  is  therefore  transformed  into  a  sinusoidal 
curve  of  180  points  in  the  accumulator,  where  each(Rp) 
element  of  the  curve  describes  possible  lines  through 
that  image  point.  If  points  in  image  space  are 
members  of  the  same  line,  curves  formed  by  those 
points  will  cross  in  parameter  space  at  the  (£.p)elerent 
which  describes  the  line  (Fig.  lb).  That  element 
will  have  a  higher  value  than  surrounding  elements  In 
the  accumulator  array  (Duda  and  Hart,  1972).  The 
accumulator  array  can  be  thought  of  as  a  2-d 
histogram  of  the  frequency  with  which  points  ocrur  on 
a  given  line.  Feaks  in  the  accumulator  above  a  noise 
threshold  correspond  to  lines  in  the  image. 


IMAGE  SPACE 
(0.0)  (n.O) 


PARAMETER  SPACE 
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Fig.  1.  The  Hough  transform  maps  points  from  image 
space  (la)  into  parameter  space  (lb).  text  f«r 
explanation. 
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3.  APPLICATION  TO  LEAD  FEATURES 


3.1  Tost  Images 

Figure  2  shows  the  AVHRR  image  of  the 
central  Arctic  from  which  test  images  1-4  were 
selected.  Image  resolution  varies  with  distance  from 
nadir.  For  simplicity,  ve  will  assume  a  constant  1km 
per  pixel.  Leads  are  varmer  than  the  surrounding 
ice,  and  therefore  appear  ns  dark  lines  in  infrared 
imagery.  Binary  images  were  created  from  Images  1-4 
using  an  interactive  detrending  and  thresholding 
procedure  (Fig.  3b,  5a).  This  preprocessing  step  Is 
critical,  aS  the  Hough  transform  will  detect  only 
those  lines  which  appear  in  the  binary  Image.  Future 
work  will  seek  to  make  this  step  less  subjective  and 
to  preserve  the  distinction  between  new,  dark  leads 
and  old,  light  leads  which  have  a  thicker  ice  cover, 
rot tunntely,  the  Hough  transform  is  insensitive  to 
breaks  in  lead  lines  which  arc  created  in  the 
thresholding  process. 


3.2  Transform  Procedure 

Fig.  A a  displays  the  accumulator  for  the  image  of 
Fig.  3l>  in  image  form.  The  resolution  of  each  pixel 
is  one  degree  in  the  0  direction,  and  ?km  in  the  P 
direction.  Peaks,  or  bright  points,  in  the 
accumulator  give  the  angular  orientation  (0)  and 
noimal  distance  to  Image  origin  (f>)  for  leads  In  the 
Image.  Because  lead  pixels  are  only  approximately 
/ulliicni,  peaks  will  lie  spread  over  a  cluster  of 
points.  An  empirically  hosen  threshold  was  applied 
to  the  accumulator,  and  points  vith  values  above  the 
threshold  (shown  in  Fig.  b )  were  inverse- 1 rans formed 
to  determine  how  well  leads  in  the  binary  imago  were 
detected  (Fig.  3c).  Parareter  space  carles  no 
information  on  the  position  of  line  end  points, 
therefore  calculated  lines  are  drawn  with  infinite 
length. 

The  Inverse  transform  of  all  points  in  clusters  with 
values  above  t lie  threshold  results  in  a  fan  of  lines 
for  each  lead  in  image  space.  A  central  point  for 
each  cluster  in  Fig.  4b  was  i nverse- t rans formed  to 
give  an  average  representation  of  each  lead  (Fig. 

3d).  All  prominent  leads  were  detected,  but  several 
short  or  less  prominent  lends  were  not,  because  they 
do  not  have  collncar  points  greater  In  number  than 
the  empiric  lly  chosen  accumulator  threshold  of  51. 

The  neat -ve: t leal  line  on  the  right  in  F<g.  3d  is  a 
false  detec- ion  -  at  least  52  points  lie  on  that  line 
in  image  spree.  While  lowering  the  threshold 
increases  t1  e  probability  of  detecting  all  leads,  the 
number  of  false  detentions  rises  also.  Using  smaller 
Images  mlti  ates  this  problem  to  3omn  degree  by 
lessening  t;  e  chance  that  the  number  of  randomly 
col  inear  po  nts  will  exceed  the  number  of  points  In 
the  shortes  lead.  However,  detection  will  always  be 
biased  towa  d  longer  (and  straighter)  leads. 

In  Fig.  5,  nverse  transforms  for  images  2-4  using 
accumulator  cluster  centers  arc  shown.  Thresholds 
wetc  1  over i  until  t he  first  occurrence  of  a  false 
detection,  and  then  raised  to  prevent  false 
detection. 

4.  RESULTS 

In  order  t-  display  parameter  space  Information  os 
lead  mien  tlon  and  spacing,  histograms  of  t  e 
frequency  w  tli  which  points  occur  on  an  Image  line 
ate  plotted  fFlg.  6).  The  orientation  histogram  is 
found  tty  sh  mlng  the  elemet  In  columns  of  an 
arcumulato  such  as  that  o  Fig.  4b,  in  which  points 
below  the  roshold  have  be  n  set  to  0.  The  notmal 
length,  or  (stance  hlstogm,  I?  found  by  summing 


Fig.  2.  N0AA-9  AVHRR  image  from  25  April, 
1986,  enhanced  to  highlight  leads. 


row  elements.  The  distance  histograms  ave  peaks 
which  are  higher  and  sharper  than  those  for 
orientation,  in  part  because  orientation  may  change 
greatly  along  a  lead  while  the  normal  distance  to  'be 
origin  slays  roughly  the  same.  Both  lead  length  and 
the  ntimbot  of  leads  having  t he  same  orientation  o- 
normal  distance  from  the  origin  contribute  to  t 
histograms  In  Fig.  6.  A  smoothed  version  of  tin 
orientation  histogram  is  plotted  as  a  represents  inn 
of  average  orientation.  It  is  of  note  that  for 
Images  1  and  2,  orientation  peaks  are  separated 
about  30  degrees.  This  is  near  the  28  degree 
intersection  angle  Marko  and  Thor  son  (1977)  found 
cltarac tei  i s t i c  ot  Arctic  leads. 

Ttie  spacing  between  peaks  In  the  distance  histogi  " 
is  an  approximate  measure  of  lead  spacing.  The 
histogram  Is  treated  as  a  spatial  series  -nd  a 
Fourier  transform  performed.  Peaks  in  the  smoothed 
amplitude  spectrum  plotted  to  the  right  o'  each 
(I  is  I  anew  histogram  correspond  to  lead  spe.-ings. 

Unlike  taking  an  average  of  lead  spneings,  tills 
method  distills  spacing  information  while  preserving 
the  bimodal  natute  of  lead  spneings  such  as  those  in 
imago  1  . 


5.  CONCLUSIONS 

The  Hough  transform  simplifies  the  analysis  of  (map's 
(nr  lend  statistics.  It  yields  lerd  ctientMlon  an  i 
spacing  I nf oi mat  I  on  quickly  and  In  a  form  which  is 
easy  tn  Interpret.  As  demonstrated  here,  It  is  an 
Incomplete  method  In  that  some  pi  e-process i nR  is 
required  to  ptoduce  a  hi  tint  y  Image,  and  post- 
pt u:-i".u;tng  must  tie  done  to  extinct  average  lead  width 
and  length.  Within  these  limitations!  the  Hough 
ti  Austin  m  penults  semi  automated  extraction  of  lead 
parameters  from  digital  imagety.  Kutuie  votk  will 
attempt  to  fully  automate  t lie  ptocedute  by  setting 
the  accumulator  threshold  using  criteria  to  avoid 
false  detections  (e.g.  net ig,  (19R?)).  In  addition, 
the  slat  1st  li.il  accittacy  of  the  method  mist  he 
establ I  shed. 
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Figmn  ft.  load  orientation  and  spacing  from  Hough 
ttansfnrm  parameter  space.  Histogtams  on  left, 
smoothed  m  initiation  histograms  and  smootlied  Fourier 
transform  of  distance  histograms  on  right. 


